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Abstract 



r"| ' It is shown that, in the presence of correlations in particle emis- 

sion, the measured HBT radii are related to the correlation range 
QhI rather than to the size of the interaction volume. Only in the case 

^ ■ of weak correlations the standard interpretation may be applicable. 

The earlier discussion ^ of the short-range correlations in configu- 
ration space is generalized to include also the correlations of particle 
^ ■ momenta. 



1. Measurements of HBT correlations in multiparticle production pro- 
vide important information on the production mechanism, in particular on 
the space-time structure of the particle emission region 0. To obtain this 
information, however, it is necessary to rely on some specific theoretical in- 
terpretation of the observed phenomena. The results are model dependent: 
The physical meaning assigned to the measured quantities does depend on 
the theoretical input. 

In the standard treatment of this problem one usually starts with a 
model where particles are uncorrelated (except for Bose-Einstein correla- 
tions) and then corrects the results by including final state interactions. This 
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includes corrections for Coulomb interactions, low energy particle interac- 
tion parametrized by scattering lengths and effects of resonances jS]. In the 
present paper we discuss correlations due to strong interactions in the pro- 
duction process. Some such correlations are known to occur [4^ , some others, 
still hypothetical, may be - hopefully - uncovered by the HBT measurements 

HI- 

To simplify the presentation we consider only the two-dimensional (trans- 
verse) distributions^, taken as Gaussians to avoid complicated integrations 
which only obscure the essential points of our argument. In this case Wigner 
functions H^(pi, ...pn,; xi, ...x„) can be used instead of the more complicated 
emission functions S{pi, ...pn', Xi, The Wigner functions are real func- 

tions of momenta and positions and are in a well-defined sense the best 
quantum analog of particle density in phase-space. Therefore the parameters 
characterizing the Wigner functions can be interpreted^ as the parameters 
characterizing the space distribution of sources and their momentum spectra 



The density matrix in momentum space is related to the Wigner function 
by the formula: 

p(pi,...,p„;p;,...,p;) = 
= J (fxi...(fxn exp [i (QiXi + ... + Q„x„)] W{Ki, K„; Xi, ...x„) (1) 

where = (pi + p-)/2 and Qi = Pi - p-. 

It follows that the momentum distribution of particles can be expressed 

as 

fio(Pl, •••,Pn) =p(pi,...,Pn;Pl,...,Pn) = 

= j ci^a;i...(i^x„iy(pi, ...,p„;xi, ...x„) (2) 
Similarly, for the coordinate distribution we have 

i7o(xi, x„) = p(xi, x„; Xi, x„) = 

= j rfV---C^^PnVr(pi, ...,p„;xi, ...x„) (3) 

"'^i.e. distributions integrated over some interval of the longitudinal variables. 

^Given all the caveats related to the fact that we are dealing with quantum phenomena 

[IE]. 
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For the momentum distribution of identical bosons we have to symmetrize 
the production amphtudes. This modifies the momentum distribution (see, 
e.g., [H]) into 

n{pi, p„) = ^ ^ p{pi^, pi^; Pi, , p,,J (4) 
^* p pf 

where the sum runs over all permutations P and P' of (ii, ...in) and (i^, ...i'n)-^ 
This is the key formula which explains the main interest in the HBT mea- 
surements: the distribution of identical particles opens a window to the non- 
diagonal elements of the density matrix and thus also to the Wigner function. 
It is also clear, however, that this information is not sufficient to obtain full 
information about the distribution of sources. Thus further theoretical input 
is needed. 

The purpose of the present paper is to discuss the physical meaning of 
the measured two-particle HBT parameters in terms of the characteristics 
of the momentum and coordinate distribution of the sources as described 
by the Wigner function. The well-known case of uncorrelated emission (for 
recent reviews, see e..g. 0) is summarized briefiy in the next section. The 
emission of particles correlated in pairs is described in Section 3. In Section 
4 a more realistic situation, when only a fraction of the particles is emitted 
in pairs while others remain uncorrelated, is considered. The experimental 
consequences are discussed in Sections 5 and 6. Our conclusions are listed in 
the last section. 

2. The assumption of uncorrelated production means that the Wigner 
function factorizes into a product of single particle Wigner functions. Of 
course this factorization is then satisfied also for the unsymmetrized density 
matrix. 

To illustrate the consequences of this Ansatz and to fix our notation, con- 
sider a single particle Wigner function in the most general Gaussian form^'^ 



iy(p,x) 



exp 



2A^ 



(x - r„p/A„) 
2(i?S - rl) 



(5) 



•^For fermions there is an extra minus sign when P andP' are odd with respect to each 
other. 

"^As aheady mentioned in the Introduction, aU vectors are two-dimensional. 



^This model is sometimes referred to as the Zajc model P|. 
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One sees that the parameter r„ is responsible for momentum-position corre- 
lation. From (0), using Q and Q, we derive for single particle distributions 



^^o(p) = Jd^xW{p,^) = ^^exp 
Qoi^) = J rfW(p,x) = ^^^exp 



2A? 



X 



2m 



(6) 



One sees that the parameter A„ describes the width of the distribution 
in momentum space whereas Ru determines the size of the system in config- 
uration space. 

Using (P) and ^ , we obtain the two-particle distribution for identical 
particles: 



^(Pl,P2) 

where 



exp 



P1 + P2 

2A2 



{1 + 



exp 



iPi 



V2YR\bT 



} (7) 



R]{BT = RI-i^u- T-r^ (8) 

u 

One sees that in this simple case measurements of the single particle 
distribution and pair distribution allow to determine A„ and Rhbt- One 
also sees from (jSI) that these two parameters are not sufficient to determine 
Ru, the size of the system in configuration space jTUI- To this end it is 
necessary to know the correlation between the momentum and the position 
of the emission point of a particle, as expressed by the parameter ru- 

3. The most general Gaussian two-particle Wigner function, symmetric 
with respect to simultaneous exchange of the particle momenta and positions, 
can be written as 



iyc(pi,P2;xi,X2) 



167r4A^A^(i?i - rl){Rl - rl] 



exp 



PI 

A2_ 



exp 



(x+ -r+p+/A+) 
Rl - rl 



[x_ — r_p_/A, 
Rl-rl 



(9) 
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where p± = (pi ± P2)/2 and xj. = (xi ± X2)/2. Note that if 



A_ = A, 



R-L- = R- 



(10) 



the Wigner function factorizes and the problem reduces to the one discussed 
in the previous section. 

One sees from Q that r± are responsible for the correlations between po- 
sitions and momenta. To see the physical meaning of the other 4 parameters 
we calculate the distribution of momenta 



^o(Pi,P2) 
and positions 



47r2A^A2 



exp 



Pi + Pj 
2AI 



IPi - P2j 

2cj2 



(11) 



fio(Xi,X2) 



where 



1 



An^RlRl 



exp 



2Rl 



2A' 



2AV 



1 

e 



1 



(Xi - X2j 
2^2 



1 



2Ri 2Ri 



(12) 



(13) 



From this we see that A^ describes the momentum distribution, whereas 
C(j2 describes the correlations between the momenta in the pair. Similarly, 
R^ describes the distribution of the particle positions while describes 
correlations between the positions of particles in the pair. Note that 
and ^2 are not necessarily positive. Note also that correlations do indeed 
disappear {l/u = 1/^ = 0) when condition ^TU^ is satisfied. 

Using 0, (HI) and (jlj), the two-particle density matrix is obtained: 



^(Pi, P2) = ^^o(pi, P2) (1 + exp [-(pi - P2)^i?c]) 
where Qq is given by (fTT|) and 



(14) 



Rl - 



1/4AL 



(15) 
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One sees that n(pi,p2) depends only on three parameters: A^,A^, and 
— r"^, whereas and do not have any impact on the momentum 

distribution. 

Using (fT^ we obtain 



which exphcitely shows the effect of correlations in configuration space on 
the physical interpretation of the HBT measurements. 

Note that for positive correlations (^^ > 0) R^ is always smaller than 
both ^^/2 and R^. In particular, when <C R^ we have R^ ^ ^'^/2. In this 
case the HBT measurements give only information on correlations and not 
on the size of the system in configuration space. 

One also sees that for negative correlations R^ is always greater than i?^. 

This discussion shows that correlations in configuration space can signi- 
ficatly influence the interpretation of the measured HBT parameters. Only 
if there are no correlations (1/^^ = 0), R+ and R- are identical and by this 
"accident" one can obtain information about the total volume of the reaction 

4. In the previous section we have discussed the situation when all pairs 
of the emitted particles are correlated. This is an interesting theoretical ex- 
ercise which, however, hardly corresponds to reality. The measured HBT 
correlations indicate that the data are in reasonable agreement with the as- 
sumption of uncorrected production. This suggests that to discuss practical 
consequences of our formalism it is more appropriate to consider a situation 
when correlated emission affects only a fraction of all the particles, the others 
remaining uncorrelated. 

The formalism developped in Sections 2 and 3 is well suited to cover this 
case. We write the Wigner function as a sum of two terms: One describ- 
ing the uncorrelated emission and the other responsible for the correlations. 
Following the discussion of sections 2 and 3 we write 



W{pi, ps; xi, xa) = ty„PV„(pi, Xi)iy„(p2, xg) + WcWdpi, P2; xi, X2) (17) 

where Wu{p, x) is given by (0) and VFc(Pi, P2; xi, X2) by ©• Wu is the prob- 
ability that the considered particles are uncorrelated and Wc = I — Wu is the 
probability that they were emitted as a correlated pair. 
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The density matrix is thus given by a sum of two terms, one constructed 
from Wu and the other from Wc- This gives the single particle momentum 
distribution^ 



^^o(Pi) = ^e-P?/2^"$o(pi) (18) 
271" 

where 



OAS 

$o(Pi) =w^ + ^c^f^e-^y^' (19) 

represents the modification of the single particle spectrum due to the corre- 
lated emission. Here 



^ ~ + Al ~ 2A2 ■ 

Using (fT7|) and employing (P) and (jH), the momentum distribution for 
identical particles f2(pi,p2) can now be derived and thus one can construct 
the usually measured quantity 



^^P^'P^^ = l^o(Pi)^^o(p.) ^''^ 

where f2o(pi) is the single-particle distribution in the events whith at least 
one pair of identical particles, given by ()18j) . The result is 



C(Pi, P2) = WuCuipi, P2) + WcCcipi, P2) (22) 

with 



1 _L ^-{pi-p2 fRjjBj, 

CuiPl, P2) = , , , (23) 
$o(pi)$o(p2) 



^In H18(l the corrections due to BE correlations are neglected. They are expected to be 
small at high energies. 
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and 



Cc(Pl,P2) 



g-(pi+p2)V2x+e-(Pi-P2)V2x^ 



<l>o(pi)$o(p2) 



(24) 



with 



5. The formulae ()22|) - P^ describe the HBT measurements for a general 
superposition of uncorrelated and correlated emission. They thus cover a 
wide range of possible physical situations. 

To discuss their interpretation we have to consider the possible origin of 
these two contributions. The uncorrelated emission may stem either from 
directly produced pions or from the pions emitted from uncorrelated clusters 
(resonances). The correlated emission may reflect (i) a genuine structure of 
the source ^ or (ii) the interaction between pions. The attractive interac- 
tions lead to positive correlations (^^ > 0). They are usually represented as 
clusters of pions. The repulsive interactions (which were never observed^) 
would give negative correlations (^^ < 0). 

As seen from ()22|) - P^ . for positive correlations one may expect the two 
components, Cu and Cc, to have different ranges in (pi — P2)^- The difference 
may be large, especially in heavy ion collisions. Indeed, in this case the 
range of the first one (~ l/Rfj^rp) is determined by the size of the whole 
system, whereas the range of the second one (~ 1/-Rc) is determined by 
the geometrical size of clusters (and/or of local fluctuations) and by the 
momentum distributions. 

We shall consider in detail the generic scenario when all particles are 
emitted from uncorrelated sources jTj. The single particle distribution is 
then fully determined by the distribution and decay properties of the emitting 
sources. The condition 

J d^X2(fp2Wc{jii, P2; xi, X2) = WuijPu xi) (26) 

^As already stated in Section 1, we discuss here only correlations due to strong inter- 
actions in the production process. 
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implies 



2Al = ^l + Al- 2Rl = Rl + Rl- 2r„A„ = r_ A_ + r+A+ (27) 

and, naturally, $o(p) = 1- 

A special case of this scenario (particle emission from independent gran- 
ules) was discussed in P where it was furthermore assumed that (i) the 
distribution of sources is momentum-independent (1/A^ = 0) and (ii) the 
momentum dependence in source decay may be neglected with respect to 
dependence on difference of momenta (1/A^ <C -R^' ^hbt)- Under these 
conditions^ the expression for the correlation function considerably simpli- 
fies 

a4 

C(pi, P2) = 1 + WuC-^P'-P^^'^HBT + yj^^-^e-^P'-P^^'^' (28) 

/_\_|_ Z_A _ 

where Wc = 1/n and n is the total number of sources. 

One sees clearly the two-component structure of the correlation function^. 
As pointed out in P, the observation of the second term may serve as an 
indication of the clustering and/or of the granular structure of the emission 
region in heavy ion collisions. The size of the granules (clusters) may be read 
off from the range of the second component. 

The simple formula ()28|) illustrates very well the basic physics of the 
problem. As seen from our general expression (j23), however, the actual shape 
of the second component may be significantly influenced by the momentum 
dependence of the emitting sources. It is true that 1/A^ and 1/A^, being 
of the order of 1 fermi^ or less, are small as compared to R^bt which (in 
heavy ion collisions) is of the order of (several fermi)^. They may well be 
comparable, however, with R^ which need not be much larger than 1 fermi^. 
Thus neglecting the momentum dependence of the emitting sources P may 
be a too drastic simplification. 

Moreover, even in absence of the correlations in configuration space (i.e., 
for i?+ = i?_;r_|. = r_ = 0) the two component structure of the correlation 
function persists. Indeed, we obtain from 

^They are too restrictive: to obtain (|28|) it is enough to assume A+ — A+ — A„, i.e., 
no correlations in momentum space. 

sum of two Gaussians in the two-particle correlation function was also considered 
for another reason in 
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C(pi, P2) =wu+ [w^ + ^c^^e-(P^+P^)'/^^+^ g-(p,-P2)^i?.W + 



L\\L\_ 



The two-component structure is recovered but now the momentum corre- 
lations and not the correlations in configuration space are responsible for 
it. 

We conclude that, although the two-component structure of the HBT 
measurements seems a robust consequence of the correlated emission, the 
physical meaning of the measured parameters is by no means unique. Thus 
we feel that in the analysis of actual experiments our general approach, sum- 
marized in the formulae ()22 |) - P^ . may be needed to account for the observa- 
tions and to give the correct physical meaning to the measured parameters. 

6. Several comments are in order. 

(i) One may note that, since for positive correlations one naturally expects 
> A^, (jZZI) implies that X- > and x+ < 0. This means that (c.f. 

(j211)) increases with increasing momentum of the pair. This effect may turn 
out helpful for identification of the second component^^. 

(ii) It is worth to remember that there are several reasons why the condi- 
tions (Hlj), relating the correlated and the uncorrelated distributions, may be 
violated (also the probability Wc of correlated emission need not be equal to 
l/n). First, not all particles are emitted in clusters, some of them are pro- 
duced directly. Second, most of the clusters observed in hadronic collisions 
are characterized by fairly small multiplicity (about three particles on the 
average) and rather small charge Therefore only a small fraction of all 
clusters emit two identical charged pions and there is no obvious reason why 
they should have the same properties as an average cluster. Thus although 
one may hope that the discussion of the previous section describes correctly 
the basic physics of the problem, the quantitative analysis may require the 
more fiexible approach. 

(iii) Finally, let us comment on the possibility of negative correlations, 
i.e. repulsive interaction (^^ < 0, tu^ < 0). In this case the cluster picture 
is not applicable. From (fT^ we deduce R- > -R+. Since _R+ is expected to 

^°This conclusion relies heavily on the condition H27() and thus needs not be generally 
valid. 
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be close to Ru, we conclude that Rc > Rhbt, i-e., the range of the second 
component is shorter than that of the first one. Thus an observation of an 
abnormally narrow peak in the distribution of (pi — P2)^ may be an indication 
of repulsive interactions in the system. It would be interesting to analyze the 
data keeping this perspective in mind^^. 

7. In conclusion, we have analyzed the effects of interparticle corelations 
in particle emission on the measurements of quantum interference. It has 
been shown that the physical interpretation of the measured parameters is 
significantly influenced by the presence of such correlations. In particular, for 
strongly correlated systems the measured range of the HBT effect is related 
to the correlation range rather than to the size of the interaction volume. 
Only in the case of weak correlations the standard interpretation may be 
applicable. The short-range positive correlations in configuration space were 
discussed in detail. The analysis given in (Ij was generalized. A possibility 
to uncover negative interparticle correlations, if any, was pointed out. 
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